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INTRODUCTION 


Waves  breaking  on  a  beach  cause  sediment  to  be  transported  both 
parallel  to  (longshore)  and  perpendicular  to  (on-offshore)  a  beach. 
Recognizing  the  importance  of  these  processes  to  Navy  operations,  the 
Independent  Research  Program  of  the  Naval  Civil  Engineering  Laboratory 
(NCEL)  has  sponsored  research  to  develop  improved  sediment  transport 
models  for  the  surfzone.  The  present  study  deals  only  with  on-offshore 
transport . 

Although  on-offshore  and  longshore  sediment  transports  are  mani¬ 
festations  of  the  same  process,  past  investigations  have  generally 
treated  them  separately  for  reasons  of  simplicity.  Fortunately,  this 
separate  treatment  has  been  relatively  successful  for  the  following  two 
reasons : 

1.  Wave  refraction  generally  causes  the  waves  to  have  near  normal 
incidence  to  the  beach  at  breaking. 

2.  The  mean  longshore  and  on-offshore  currents  are  generally  much 
weaker  than  the  oscillatory  velocity  magnitude. 

As  a  result,  the  longshore  transport  has  been  found  to  be  relatively 
well  modeled  (e.g.,  Inman  and  Bagnold,  1963;  Komar,  1971)  as  the  product 
of  an  oscillatory  velocity- induced  sediment  load  and  a  transport  veloc¬ 
ity  proportional  to  the  longshore  current.  Similarly,  the  on-offshore 
transport  may  be  modeled  as  a  balance  between  gravity,  the  asymnetry  of 
the  oscillatory  velocity  distribution,  and  the  on-offshore  steady  current 
(e.g.,  Inman  and  Frautschy,  1966;  Bowen,  1980;  Bailard,  1981). 

Although  the  above  processes  have  been  qualitatively  understood  for 
some  time,  relatively  greater  progress  has  been  made  in  quantitatively 
predicting  the  longshore  transport  as  opposed  to  predicting  the  on- 
offshore  transport.  Several  aspects  of  the  latter  process  tend  to  make 
modeling  efforts  more  difficult.  These  include: 
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1.  A  breaking  wave/bore  propagation  model  suitable  for  predicting 
oscillatory  velocity  asymmetries  and  mean  on-offshore  currents  inside 
the  surf zone  has  been  lacking. 

2.  The  effect  of  the  downs lope  component  of  the  sediment  load  on 
the  on-offshore  sediment  transport  has  generally  not  been  adequately 
represented  in  the  sediment  transport  model  used. 

3.  The  dynamic  near-equilibrium  of  on-offshore  transport  on  a 
stable  beach  tends  to  magnify  the  importance  of  second  order  effects 
such  as  the  vertical  velocity  structure  or  the  threshold  criteria  for 
sediment  movement. 

Because  of  these  difficulties,  a  number  of  past  models  have  sought 
to  correlate  on-offshore  sediment  transport  with  average  incident  wave 
conditions,  without  considering  many  of  the  details  of  the  fluid-sediment 
motions.  Some  examples  of  this  type  of  approach  are  discussed  in  Seymour 
and  King  (1982)  and  include:  the  wave  steepness  models  by  Dean  (1973) 
and  Hattori  and  Kawamoto  (1981);  the  wave  height  models  by  Saville 
(1957),  Aubrey  (1978)  and  Short  (1978);  and  the  wave  power  model  by 
Short  (1978).  A  few  models  have  included  some  degree  of  detail  of  the 
fluid  sediment  motions.  Examples  of  models  which  consider  wave  velocity 
asymmetry,  downslope  sediment  loads,  and  mean  on-offshore  currents 
include  those  by  Inman  and  Frautschy  (1966),  Bowen  (1980)  and  Bailard 
(1981).  These  three  models  are  based  on  adaptations  of  Bagnold's  (1963, 
1966)  sediment  transport  model  for  streams.  Bailard's  (1981)  model  is 
the  most  complete  in  that  nonnormal  wave  incidence  is  permitted  as  is 
the  presence  of  longshore  currents . 

A  common  aspect  of  Bowen's  (1980)  and  Bailard's  (1981)  energetics- 
based  models  is  the  importance  of  several  surfzone  velocity  moments  in 
determining  the  direction  and  magnitude  of  the  on-offshore  sediment 
transport.  These  moments  are  defined  in  terms  of  idealized  monochromatic 
waves.  However,  they  can  be  extended  to  spectral  wave  inputs  using  the 
techniques  described  in  Guza  and  Thornton  (in  review).  Using  Stokes' 
second-order  wave  solution  and  Longuet-Higgins '  (1953)  bottom  streaming 
model  to  estimate  the  wave  velocity  asymmetry  and  mean  on-offshore 
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current,  respectively,  Bowen  (1980)  and  Bailard  (1981)  were  able  to 
qualitatively  describe  the  equilibrium  beach  profile  as  a  function  of 
the  incident  wave  amplitude  a;  the  wave  frequency  o;  and  sediment  fall 
velocity  W  (Figure  1).  Small  amplitude,  long  period  waves  were  found  to 
produce  a  steep  beach  while  large  amplitude,  short  period  waves  were 
found  to  produce  a  flat  beach.  Large  diameter  sand  grains  with  high 
fall  velocities  were  also  observed  to  produce  steeper  beaches  than 
smaller  diameter  sand  grains  with  lower  fall  velocities.  These  results 
qualitatively  confirm  observed  beach  behavior,  but  neither  the  nonlinear 
wave  solution  nor  the  mean  current  solutions  are  valid  inside  the  surfzone. 
Unfortunately  other  wave  shoaling  models  are  equally  invalid  inside  the 
surfzone  due  to  the  dissipation  of  energy  and  the  generation  of  turbulence 
associated  with  the  breaking  wave. 

Because  of  the  inability  of  existing  wave  shoaling  models  to  accu¬ 
rately  describe  wave  velocity  asymmetry  and  mean  on-offshore  currents 
inside  the  surfzone,  little  is  known  about  these  qualities.  Two  recent 
studies  have  described  limited  field  measurements  of  these  quantities 
(Huntley  and  Bowen,  1975;  Guza  and  Thornton,  in  review).  These  studies 
suggest  that  the  magnitudes  of  these  quantities  may  be  a  function  of  the 
incident  wave  conditions,  the  beach  slope,  and  the  local  water  depth. 

Because  of  a  lack  of  good  field  data,  very  few  comparative  evalua¬ 
tions  have  been  made  of  existing  on-offshore  sediment  transport  models. 
Recently,  however,  a  series  of  large-scale  field  experiments  have  been 
conducted  as  part  of  the  Nearshore  Sediment  Transport  Study  (NSTS).  The 
first  of  these  experiments  was  conducted  at  Torrey  Pines  Beach,  Calif., 
during  November  1978.  During  this  time,  simultaneous  measurements  of 
deepwater  wave  characteristics,  surfzone  nearbottoro  velocity  distribu¬ 
tions,  and  beach  profile  changes  were  measured.  Details  of  the  experi¬ 
ment  are  described  in  Gable  (1979).  Seymour  and  King  (1982)  utilized 
this  data  set  to  evaluate  a  variety  of  on-offshore  sediment  transport 
models.  Although  none  of  the  models  showed  great  ability  to  predict  the 
measured  beach  volume  changes,  the  simple  models  utilizing  incident  wave 
height,  wave  power,  or  wave  steepness  could  account  for  between  25%  to 
35%  of  the  total  variance  in  the  beach  volume  measurements.  The  only 
energetics-based  model  evaluated  was  a  bedload  model  developed  by  Bailard 
and  Inman  (1981)  which  performed  rather  poorly  by  comparison. 
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For  two  reasons  the  poor  performance  of  Ballard  and  Inman's  bedload 
model  is  not  surprising.  First,  studies  by  Bowen  (1980)  and  Bailard 
(1981)  have  suggested  that  most  of  the  transport  in  the  surfzone  is  by 
nearbottom  suspension.  Second,  the  velocity  moments  controlling  on- 
offshore  sediment  movements  in  Bailard  and  Inman's  (1981)  bedload  model 
are  very  sensitive  to  infrequent  data  errors  that  are  present  on  the 
tapes.  Without  a  careful  prescreening  of  the  data,  the  estimated  on- 
offshore  sediment  transports  must  be  viewed  as  suspect. 

The  success  of  Bowen's  (1980)  and  Ballard's  (1981)  total  load  surf- 
zone  models  in  qualitatively  describing  the  variation  of  the  equilibrium 
beach  slope  with  wave  and  sediment  characteristics  suggested  that  Bagnold's 
energetics  approach  may  have  some  ability  to  quantitatively  predict 
on-offshore  sediment  transport.  Consequently,  a  study  was  initiated  to 
explore  the  ability  of  Ballard's  (1981)  total  load  sediment  transport 
model  to  predict  daily  beach  volume  changes  using  data  from  the  NSTS 
experiment  at  Torrey  Pines  Beach,  Calif.,  in  November  1978.  The  objec¬ 
tives  of  the  study  were  threefold.  First,  the  model  was  tested  in  its 
ability  to  predict  daily  beach  volume  changes  using  relatively  short  (1- 
to  2-hour)  surfzone  current  meter  records.  Second,  a  relationship  was 
sought  between  average  values  of  surfzone  velocity  moments  and  incident 
wave  characteristics.  Lastly,  these  average  surfzone  velocity  moments 
were  used  to  predict  on-offshore  sediment  transport  as  a  function  of 
significant  wave  height. 

SEDIMENT  TRANSPORT  MODEL 

Ballard's  (1981)  total  load  sediment  transport  model  is  based  on  an 
adaptation  of  Bagnold's  (1963,  1966)  energetics-based  total  load  sediment 
transport  model  for  streams.  The  latter  is  generalized  for  time-varying 
flow  over  an  arbitrarily  sloping  bottom,  resulting  in 


where 


=  instantaneous  sediment  transport  rate  vector 
p  =  density  of  water 

=  drag  coefficient  of  the  bed 
Eg  =  the  bedload  efficiency  factor 
<|)  =  internal  angle  of  friction  of  the  sediment 
tan  p  =  the  bed  slope 

Eg  =  suspended  load  efficiency  factor 
W  =  is  the  fall  velocity  of  the  sediment 
u^  =  instantaneous  nearbottom  fluid  velocity  vector 
i  =  unit  vector  directed  upslope 
<  >  =  time-average 

Note  that  for  both  the  bedload  (first  bracketed  quantity)  and  the  sus¬ 
pended  load  (second  bracketed  quantity) ,  the  transport  consists  of  a 
primary  component  directed  parallel  to  the  instantaneous  fluid  velocity 
vector  and  a  secondary  component  directed  downslope.  The  latter  is 
associated  with  the  downslope  gravity  component  of  the  sediment  load. 

Figure  2  depicts  a  plane  contour  beach  with  the  x-axis  directed 
shoreward  and  normal  to  the  beach  and  the  y-axis  directed  parallel  to 
the  beach.  The  slope  of  the  beach  is  tan  p.  Considering  only  the 
on-offshore  sediment  transport.  Equation  1  becomes 


Equation  2  is  assumed  to  be  valid  for  any  nearbottom  velocity  field  and 
can  be  directly  used  with  data  from  a  two-axis  (x  and  y)  current  meter 
to  predict  the  time-averaged  on-offshore  transport  rate  at  the  current 
meter  location. 


For  modeling,  it  is  convenient  to  use  a  monochromatic  wave  repre¬ 
sentation  for  the  nearbottom  water  velocity  field.  Bailard  (1981) 
assumed  a  velocity  field  composed  of  an  oscillatory  velocity  component  u 
oriented  at  an  angle  a  to  the  x-ax.i.s  and  a  steady  velocity  component  u 
oriented  at  an  angle  6  to  the  x-axis  (Figure  2).  The  total  velocity 
vector  u^  then  becomes 

=  (u  cos  a  +  u  cos  0)  i  +  (u  sin  a  +  u  sin  0)  j  (3) 


In  addition,  the  oscillatory  velocity  component  is  assumed  to  be  asym¬ 
metrical  being  composed  of  a  primary  component  u^  with  frequency  a  and  a 
secondary  harmonic  u^2  with  frequency  2o  so  that 

u  ~  u  cos  at  +  u  -  cos  2at  +  ...  (4) 

m  m2 


Substituting  Equations  3  and  4  into  Equation  2,  Bailard  (1981)  obtained 
the  idealized  on-offshore  transport  equation 
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For  weak  mean  currents,  Snell's  law  and  the  spilling  wave  hypothesis 
are  sufficient  for  obtaining  estimates  of  a,  u^,  u*,  and  u*  throughout 
the  surfzone  as  a  function  of  incident  wave  conditions. t  Moreover,  a 
longshore  current  model  (e.g.,  Longuet-Higgins ,  1970)  is  sufficient  for 
estimating  the  magnitude  of  6^.  Unfortunately,  the  two  skewness  para¬ 
meters  *1',  and  t-,  as  well  as  the  normalized  mean  onshore  current  6  , 
12  u 

cannot  be  estimated  from  present  surfzone  wave  shoaling  models.  Well 
outside  the  surfzone,  where  a  rippled  bottom  generally  invalidates  the 
1  .’esent  sediment  transport  model,  and  6^  may  be  estimated  from 

Stokes'  second-order  wave  solution  and  Longuet-Higgins’  (1953)  bottom¬ 
streaming  solution.  Assuming  negligible  longshore  currents,  weak  on- 
offshore  currents  and  normal  incidence,  Bowen  (1980)  and  Bailard  (1981) 
used  these  solutions  to  obtain  the  previously  discussed  solution  for  the 
equilibrium  beach  slope. 

Neither  Stokes'  second-order  wave  solution  nor  Longuet-Higgins' 
(1953)  bottom-streaming  solution,  however,  are  valid  inside  the  surfzone. 
Other  wave  shoaling  models  are  invalid  there  as  well.  One  of  the  objec¬ 
tives  of  this  study  was  to  examine  measured  surfzone  values  of  “Vj ,  H'^, 
and  6^  using  data  from  the  NSTS  Torrey  Pines  experiment  (Gable,  1979). 

tNote,  however,  that  Guza  and  Thornton  (in  review)  found  that  the 

commonly  assumed  spilling  wave  hypothesis  did  not  lead  to  an 

accurate  estimate  of  u  . 

m 
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In  this  respect,  the  present  study  was  a  continuation  of  the  study  by 

Guza  and  Thornton  (in  review),  who  showed  that  the  quantities  defined  in 

Equation  5  are  meaningful  only  for  monochromatic  waves  incident  from  a 

single  direction  a.  On  an  actual  beach  the  incident  waves  compose  a 

spectrum  with  varying  energy  content  at  different  frequencies  and  wave 

angles.  As  a  result,  Guza  .d  Thornton  (in  review)  defined  equivalent 

quantities  for  *?, ,  f „ ,  u  ,  u*  and  u*,  which  may  be  estimated  from  the 
1  z  m  j  b 

measured  current  meter  data.  In  addition,  a  series  of  six  equivalent 
wave  angles.  Of,  to  a,,  are  defined  for  the  different  quantities  in 
Equation  5.  Equation  5  then  becomes  (Guza  ai.d  Thornton,  in  review) 
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~  2  ~2^~2  ^  -2  -2^-2 

u^  =  u  +  V  and  u.j.  =  u  +  v 
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In  addition,  the  equivalent  angles  ,  a^t  calculated 

as  follows 


tan  Oj 


~-3  ~2  ~ 

<v  *  u  v> 

~3  ^  ~  ~2^ 

<U  +  U  V  > 


(19) 

(20) 

(21) 

(22) 


Guza  and  Thornton  (in  review)  analyzed  2  days  of  NSTS  Torrey  Pines  data 
and  presented  estimated  values  for  4*^,  “3*  “5>  ®j*  ®2*  ®3’ 

and  or^  as  a  function  of  depth.  They  found  that  for  the  conditions 
studied  6  «1,  6  «1,  and  a-  and  a.  were  small.  In  addition,  the 

present  study  has  shown  generally  that  cos  "  1  and  cos  ~  ±1.  As  a 
result  for  conditions  such  as  those  at  Torrey  Pines  Beach,  Equation  13 
can  be  simplified  to 


where  and  are  assumed  to  be  equal  to  4*^  cos  and  4*2  cos 
respectively,  and  can  thus  take  on  negative  (offshore-directed)  values. 

Using  Equation  23,  estimating  the  on-offshore  sediment  transport 
rate  becomes  a  relatively  simple  matter  of  estimating  u^,  4*^,  4*2,  6^, 
u^,  and  u*.  One  of  the  objectives  of  the  present  study  was  to  determine 
whether  these  parameters  can  be  estimated  from  average  gross  character¬ 
istics  of  the  incident  waves . 
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DATA  ANALYSIS 


Data  Set  Description 

In  November  1978,  a  month-long  field  experiment  was  conducted  at 
Torrey  Pines  Beach  as  part  of  the  Nearshore  Sediment  Transport  Study 
(NSTS) .  Torrey  Pines  Beach,  Calif.,  is  a  plane-contoured  beach  with  a 
concave  profile.  The  slope  of  the  beach  face  is  approximately  0.05, 
decreasing  to  about  0.02  within  the  surf zone  (Figure  3).  The  sand  on 
the  beach  is  moderately  well  sorted  with  a  mean  diameter  of  0.17  mni. 

The  beach  exhibits  no  bar-trough  features.  The  experiment  consisted  of 
the  simultaneous  measurement  of  incident  wave  conditions,  nearbottom 
water  velocities  distributions,  sand  tracer  movements,  and  beach  profile 
changes.  The  incident  wave  climate  was  measured  in  10  meters  of  water 
with  a  linear  array  of  pressure  sensors.  Nearbottom  water  velocities 
were  measured  using  dual-axis  (x  and  y)  electromagnetic  current  meters, 
while  the  beach  profiles  were  measured  with  a  rod  and  transit  onshore 
and  a  fathometer  offshore. 

The  current  meters  were  placed  in  a  cross-shaped  pattern  within  the 
nearshore  area.  Referring  to  Figure  4,  7  current  meters  formed  a  line 
perpendicular  to  the  beach  and  ranging  in  depth  from  0.25  meters  to 
6  meters  relative  to  mean  sea  level  (MSL).  Ten  other  current  meters 
formed  a  line  parallel  to  the  beach  at  a  depth  of  about  1  meter  relative 
to  MSL.  The  beach  profiles  were  measured  at  five  ranges  along  the 
beach.  Further  details  of  the  measurements  and  the  experiment  can  be 
found  in  Gable  (1979). 

Model  Evaluation 


One  of  the  objectives  of  the  present  study  was  to  evaluate  the 
ability  of  Equation  2  to  predict  daily  beach  volume  changes  using  NSTS 
data.  The  approach  used  to  evaluate  Equation  2  was  the  same  as  that 
used  by  Seymour  and  King  (1982).  The  daily  beach  volume  changes  were 
computed  by  integrating  the  beach  profile  changes  across  100  meters  of 
beach  face.  The  seaward  extremity  of  the  integration  interval  roughly 
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coincided  with  the  location  of  the  shore-parallel  current  meters  so  that 
a  simple  box  model  analysis  could  be  used  to  estimate  the  beach  volume 
changes.  In  the  present  study  the  flux  of  sand  into  or  out  of  the  box 
was  predicted  by  Equation  2  using  the  measured  current  meter  data  from 
each  shore  parallel  current  meter.  The  measured  volume  changes  used  in 
the  present  study  are  those  tabulated  by  Seymour  and  King  (1982). 

Equation  2  has  3  free  parameters:  the  bed  friction  coefficient  c^, 
the  bedload  efficiency  factor  e„,  and  the  suspended  load  efficiency 

D 

factor  Cg.  In  the  present  study,  a  value  of  0.005  was  selected  for  c^, 
based  on  an  analysis  of  longshore  current  data  at  Silver  Strand  Beach, 
Calif.,  (Bailard,  1981).  This  beach  has  a  bed  slope  similar  to  that  at 
Torrey  Pines  Beach  (0.034)  and  is  located  15  miles  to  the  south.  On  the 
other  hand,  Thornton  and  Guza  (1982)  examined  longshore  currents  at 
Torrey  Pines  Beach  and  concluded  chat  c^  =  0.01  ±  0.01.  As  evidenced  by 
the  large  uncertainty  interval,  Thornton  and  Guza  concluded  that  the 
Torrey  Pines  data  set  was  not  particularly  good  for  estimating  the  size 
of  c^. 

The  bedload  and  suspended  load  efficiency  factors  and  are  the 
remaining  two  free  parameters  in  Equation  2.  Rewriting  Equation  2  to 
isolate  these  factors,  we  obtain 


‘S"  =  'b  'a  *  ‘s  'b  - 's  'c 


(24) 


where 


A  tan  0 


,  2  _  tan  P  -  3. 

<  u,_  u^*i>  -  - - r  <  u^  > 

t  t  tan  ^  t 


P  ‘^f  +  1 3  -»  ^ 

Ib  =  —  <  “tl  “t*^^ 


(25) 


(26) 


P  =4 


-»  5 

tan  P  <  u^  > 


(27) 


Equation  24  expresses  the  on-offshore  transport  rate  in  terms  of  immersed 
weight  transport.  The  volumetric  transport  rate  may  be  related  to 
the  immersed  weight  transport  rate  <i^^  by  the  following  equation 
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<i  > 

X 


(28) 


(P  -P)  g  N 


where  =  density  of  the  sand  grains 

g  =  gravity 

N  =  "at  rest"  volume  concentration  of  sediment  assumed  here  to 
°  be  0.6 


The  procedure  used  to  estimate  the  bedload  and  suspended  load  effi¬ 
ciencies  Cq  and  c  was  as  follows.  First,  daily  estimates  of  I.,  and 
were  obtained  from  64-minute  records  for  each  of  the  10  shore-parallel 
current  meters.  These  estimates  were  averaged  together  to  form  a  single 
daily  estimate  that  was  assumed  to  be  representative  of  the  general 
experiment  area.  Table  1  contains  a  summary  of  the  estimated  values  for 
Ig,  and  Ij..  Next,  a  correlation  analysis  was  conducted  to  determine 
the  lag  time  between  the  predicted  volume  changes  and  the  measured 
changes.  To  do  this,  previously  estimated  values  of  =  0.21  and 
Cg  =  0.025  obtained  from  longshore  transport  data  (Bailard,  1981)  were 
used  with  Equation  24.  Follovi.og  the  rationale  of  Seymour  and  King 
(1982),  only  lag  times  of  zero  and  1  day  were  considered.  The  maximum 
correlation  (R*  =  0.19)  occurred  with  a  lag  time  of  1  day,  as  was  found 
by  Seymour  and  King  (1982)  for  a  number  of  other  models  using  the  same 
data  set.  The  lag  time  result  was  relatively  insensitive  to  different 
values  of  e  and  c  ,  although  varied  as  would  be  expected. 

D  O 

The  relatively  small  value  of  R^  (0.19)  means  that  only  19%  of  the 
observed  variance  in  the  measured  beach  volume  changes  is  accounted  for 
by  the  model.  By  comparison,  using  the  same  data  set,  Seymour  and  King 
(1982)  found  R^  values  as  high  as  0.35  for  a  simple  model  relating  the 
wave  steepness  to  on-offshore  sediment  movements.  Possible  reasons  for 
the  relatively  poor  performance  of  the  present  model  are  discussed 
below. 

A  nonlinear  least-squares  estimation  procedure  (Draper  and  Smith, 
1966),  was  used  to  estimate  Cg  and  Cg  from  the  data  in  Table  1.  The 
estimation  procedure  required  constructing  a  contour  map  of  the  mean 
square  error  S,  defined  as 
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1  ^  2 

S  (e_  £_)  =  V'  (V  -  V  .  )  (29 

B  S  n-2  2^  meas  pred 

i=l 

where  n  =  number  of  data  pairs  (17) 

V 

meas  =  measured  beach  volume  change  with  a  1-day  lag  time 
Vpred  “  predicted  beach  volume  change 

Figure  5  shows  a  plot  of  S  versus  e_  and  e_.  The  minimum  mean  square 

D  b 

error  occurred  at  £„  =  0.10  and  £„  =  0.02.  These  estimates  are  quite 

O  b 

similar  in  size  to  those  estimated  by  Bailard  (1981)  (£„  =  0.21  and 

D 

£g  =  0.025)  based  on  longshore  transport  data.  Figure  5  also  shows  that 
95%  confidence  limits  on  £_  and  £„.  The  limits  are  much  broader  than 

D  b 

for  those  in  Bailard  (1981)  and  reflect,  in  part,  the  low  degree  of 
correlation  between  the  measured  and  predicted  beach  volume  changes.  In 
both  cases,  however,  the  predicted  values  of  £_  and  £„  fall  within  each 

D  b 

other's  areas  of  uncertainty. 

The  low  degree  of  correlation  between  measured  and  predicted  beach 
volume  changes  raises  questions  about  the  ability  of  Bailard's  (1981) 
sediment  transport  model  to  predict  on-offshore  sediment  transports.  It 
should  be  recognized,  however,  that  predicting  on-offshore  sediment 
transport  rates  is  a  severe  test  of  a  model  because  the  net  transport 
represents  a  small  difference  between  two  relatively  large  instantaneous 
on  and  offshore  sediment  transports.  Small  biases  in  either  the  onshore 
or  offshore  direction  can  significantly  influence  the  predicted  direc¬ 
tion  and  magnitude  of  the  net  transport.  In  addition,  a  sensitivity 
analysis  was  performed  whereby,  for  selected  days,  estimates  of  the 
on-offshore  transport  rate  were  calculated  for  consecutive  6A-minute 
segments  of  time.  For  a  day  with  moderate  waves  (November  4),  the 
standard  deviation  of  these  estimates  was  22%  of  the  predicted  mean  for 
four  consecutive  segments.  For  a  day  with  large  waves  (November  12), 
this  value  increased  to  46%  for  three  consecutive  segments.  In  addition, 
individual  current  meters  showed  a  much  wider  range  of  variation  (up  to 
800%)  between  consecutive  hours.  These  results  suggest  that  although 
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the  beach  may  generally  be  in  a  state  of  near  dynamic  equilibrium,  small 
changes  in  the  incident  wave  conditions  or  the  tide  may  strongly  influ¬ 
ence  short-term  local  on-offshore  movements.  As  a  result,  for  the  data 
set  studied,  short  current  meter  records  of  1  to  4  hours  may  not  be  of 
sufficient  duration  to  estimate  the  daily  beach  volume  changes.  Alter¬ 
natively,  the  model  may  simply  be  unable  to  estimate  on-offshore  sedi¬ 
ment  movements.  The  present  data  set  is  insufficient  to  resolve  this 
question. 

Velocity  Moment  Magnitudes 

The  above  discussion  concerned  the  ability  of  Equation  2  to  predict 
on-offshore  sediment  transports.  Equation  2,  however,  requires  a  full 
knowledge  of  the  surfzone  velocity  field  and  is,  therefore,  not  very 
useful  for  modeling.  Equation  23  is  a  greatly  simplified  version  oi 
Equation  2,  which  is  potentially  more  useful  for  modeling  on-offshore 
sediment  transports.  Unfortunately,  this  equation  still  contains  a 
number  of  surfzone  velocity  moments,  about  which  little  is  known.  These 
include  the  wave  velocity  skewness  parameters  and  (Equations  9  and 
10),  the  normalized  onshore  current  6^  (Equation  7),  and  the  normalized 
velocity  magnitudes  u*  and  u*  (Equations  11  and  12). 

As  previously  discussed,  only  u^,  u*,  and  u*  can  be  estimated  using 
linear  wave  theory.  For  normally  incident  monochromatic  waves  and  weak 
mean  currents,  u*  and  u*  are  equal  to  0.424  and  0.339,  respectively, 
while  for  a  Gaussian  input  with  the  same  total  variance  their  values 
become  0.562  and  1.13,  respectively  (Guza  and  Thornton,  in  review). 

Guza  and  Thornton  found  that  a  Gaussian  distribution  was  most  nearly 
correct . 

Assuming  wave  saturation  and  spilling  waves,  linear  wave  theory 
suggests  that 


“m  =  2 


where  y  =  H/h 

h  =  local  water  depth 
H  =  local  wave  height 


(30) 


Equation  30  suggests  that  u  should  decrease  from  a  maximuis  at  the 

m 

breakpoint  to  zero  at  the  beach.  In  fact,  field  measurements  show  (Guza 
and  Thornton,  in  review)  that  u  was  almost  constant  across  the  surfzone, 

ID 

due  to  the  presence  of  low  frequency  surf  beat  motions  within  the  inner 

part  of  the  surfzone.  The  remaining  parameters  in  Equation  23  (Vj,  ^ , 

and  6  )  are  zero  for  linear  waves  but  nonzero  for  nonlinear  waves, 
u 

Measurements  (Huntley  and  Bowen,  1975,  Guza  and  Thornton,  in  review) 
have  shown  that  they  are  in  fact  nonzero  under  actual  field  conditions 
and  may  vary  in  magnitude  and  sign  with  varying  incident  wave  condi¬ 
tions  . 

One  of  the  objectives  of  the  present  study  was  to  investigate  the 
magnitudes  of  these  surfzone  parameters  using  the  NSTS  Torrey  Pines  data 
set  and  determine  whether  their  magnitudes  and  distributions  across  the 
surfzone  could  be  estimated  from  incident  wave  conditions.  Values  for 
the  parameters  6^,  u*  and  u*  were  estimated  from  64-minute-long 

records  for  each  of  the  seven  shore  perpendicular  current  meters  closest 
to  shore.  Because  of  the  variability  of  these  quantities,  mean  surfzone 
values  were  obtained  by  averaging  the  results  of  all  seven  current 
meters.  Table  2  is  a  summary  of  these  average  surfzone  velocity  param¬ 
eters  and  the  incident  wave  characteristics  for  each  of  the  9  days 
investigated.  The  wave  characteristics  were  obtained  from  Guza  and 
Thornton  (1980;  in  review)  who  analyzed  pressure  sensor  records  measured 
in  10  meters  of  water.  In  addition,  Seymour  and  King  (1982)  reported 
significant  deep  water  wave  steepnesses  for  each  day.  These  too  are 
tabulated  in  Table  2. 

Linear  and  second-order  wave  theory  suggests  that  the  even  velocity 

moments  u^,  u*,  and  u*  should  be  related  to  the  significant  wave  height, 

while  the  odd  velocity  moments  H*,  ,  M*-,  and  6  should  also  be  related  to 

12  u 

the  deep  water  wave  steepness  and  the  average  beach  slope.  Recent 
studies  of  surfzone  wave  dynamics  by  Wright  et  al.  (1978,  1982),  Guza 
and  Thornton  (1982),  Bowen  (1980),  Huntley  and  Bowen  (1975)  and  others 
suggest  that  long  period  surfzone  motions  due  to  edge  waves,  surf  beat 
and  other  infra-gravity  wave  motions  may  influence  nearshore  velocity 
fields  and  the  resulting  nearshore  topography.  The  amplitudes  of  these 
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motions  are  thought  to  be  sensitive  to  a  surf  similarity  parameter  C 
introduced  by  Battjes  (1974),  Guza  and  Inman  (1975),  and  others.  The 
parameter  is  defined  as 


2  y'  2 

£  =  a  tan  P 

where  a  =  half  the  breaking  wave  height 


(31) 


For  small  values  of  e,  less  than  2.0  to  2.5,  reflective  conditions  with 
surging  breakers  are  observed  and  subharmonic  edge  wave  generation  may 
be  present.  For  larger  values  of  £,  dissipative  conditions  should 
prevail,  and  surf  beat  may  be  present. 

For  the  conditions  of  Torrey  Pines  Beach,  dissipative  conditions 
prevailed,  with  spilling  waves  and  a  wide  surfzone.  The  range  of  signi¬ 
ficant  wave  heights  varied  from  55  to  140  cm;  however  the  peak  wave 
period  varied  much  less,  ranging  from  approximately  10  to  14  seconds. 

The  variation  in  mean  beach  slope  during  the  experiment  was  also  very 
small,  so  that  the  slope  was  approximately  a  constant  0.02.  As  a  result, 
the  deep  water  wave  steepness  H„/L  and  the  surf  similarity  parameter  e 
were  primarily  a  function  of  the  significant  wave  height  Hg.  This 
suggested  that  for  the  Torrey  Pines  data  set,  a  relationship  might  be 
sought  between  all  of  the  relevant  surfzone  velocity  moments  and  the 

significant  wave  height  H  . 

s 

Figures  6  through  11  show  plots  of  mean  surfzone  values  of 
^u’  "^m’  ^*5  Lacking  more  data,  a  linear  relationship 

was  sought  between  variables;  the  linear  regression  lines  are  shown  in 
each  figure.  These  lines  are  described  by  the  following  equations: 


2 

0.303  -  0.00144  H 

1 

s 

0.603  -  0.00510  H 

2 

s 

6 

0.458  -  0.00157  H 

u 

s 

u 

51 

31.9  +  0.403  H 

m 

s 

u* 

• 

0.548  +  0.000733  H 

3 

s 

u* 

- 

1,50  +  0.00346  H 

5 

s 

where  u 

m 

is 

measured  in  cm/sec  and  H 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

in  cm. 
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Figures  6  through  11  also  suggest  that  for  the  conditions  found  at 
Torrey  Pines  Beach,  4*^ ,  and  6^  decrease  markedly  with  increasing  wave 
height.  This  is  in  direct  contrast  to  Stokes'  second-order  wave  solution 
and  Longuet-Higgins '  (1953)  bottom-streaming  solutions,  which  predict 
increasing  values  of  these  variables  for  increasing  wave  height.  The 
remaining  parameters,  u^,  u^,  and  u*,  behaved  more  as  would  be  predicted 
by  linear  wave  theory.  The  mean  orbital  velocity  magnitude  u^  was  found 
to  increase  with  increasing  wave  height,  while  u^  and  u^  were  relatively 
constant.  Typical  values  of  u*  and  u*  were  0.6  and  1.2,  respectively 
(very  close  to  the  theoretical  values  based  on  a  Gaussian  wave  distribution) . 

Velocity  Moment  Distributions 

While  the  average  surfzone  magnitudes  of  the  different  velocity 
moments  discussed  above  are  critical  in  determining  general  on-offshore 
sediment  movements,  spatial  details  of  these  movements  are  reflected  by 
the  distribution  of  the  moments  across  the  surfzone.  In  order  to  examine 
these  distributions,  it  was  first  necessary  to  normalize  the  moments  by 
appropriate  quantities.  For  the  variables  4'j ,  6^,  u^,  u^,  and  ug,  the 
mean  surfzone  value  of  each  variable  was  an  appropriate  normalizing 
quantity.  The  skewness  parameter  was  more  difficult  because  its  mean 
value  changed  sign.  As  a  result,  a  new  variable  4*2  defined  as 
follows : 


(38) 


This  variable  was  then  normalized  by  its  standard  deviation  and  defined 
as 


where  n  =  number  of  days  examined 

Figure  12  shows  a  plot  of  4*2sd  versus  significant  wave  height. 


(39) 
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Figures  13  through  18  are  plots  of  the  surfzone  distribution  of  the 

normalized  parameters  4*,,  4*',  5  ,  u  ,  u*.  and  u'^  as  a  function  of  the 

1  2  u*  m’  3  5 

normalized  surfzone  position  x^,  which  is  defined  as  the  distance  from 

shore  divided  by  the  nominal  surfzone  width.  For  the  present  study,  the 

surfzone  width  was  taken  as  the  distance  from  the  breakpoint  to  the 

still  water  shoreline.  Moreover,  the  breakpoint  depth  was  assumed  to  be 

equal  to  H  /0.78.  A  cursory  examination  of  Figures  13  to  18  reveals 

that  there  is  considerably  more  scatter  in  the  plots  of  4*2  and  6^  than 

in  the  plots  of  the  other  variables.  These  two  parameters  are  closely 

related  and  the  variability  in  4*'  is  believed  to  primarily  be  a  reflec- 

tion  of  the  variability  in  6^.  As  previously  mentioned,  the  mean  on- 

offshore  current  was  the  most  variable  surfzone  quantity  examined.  This 

was  again  the  case  in  the  surfzone  distributions. 

Considering  the  surfzone  distributions  of  4*, ,  u  ,  u*.  and  u* 

1  m  3  5 

(Figures  13,  and  16  through  18),  each  variable  appears  to  have  a  pre¬ 
dictable  distribution.  Whereas  4*^,  u*,  and  u*  increase  with  increasing 

distance  from  shore,  u  is  seen  to  decrease  with  increasing  distance 

m  .s 

from  shore.  The  data  for  4*2  and  6^  are  so  scattered  that  no  definite 
trends  are  evident. 


ON-OFFSHORE  TRANSPORT  SIMULATION 

The  simplified  on-offshore  sediment  transport  model  Equation  23  was 
found  to  contain  several  surfzone  velocity  moments  about  which  little 
was  known.  These  moments  included  the  wave  velocity  skewness  parameters 
and  4*2,  the  mean  on-offshore  current  6^,  the  orbital  velocity  magnitude 
u^,  and  the  normalized  total  velocity  magnitudes  u^  and  ug.  For  the 
wave  and  beach  conditions  present  during  the  NSTS  Torrey  Pines  experiment, 
average  surfzone  values  for  those  parameters  were  found  to  be  linearly 
related  to  the  significant  wave  height  through  Equations  32  to  37. 
Combining  Equation  23  with  Equations  32  to  37,  the  average  surfzone 
on-offshore  transport  rate  and  direction  can  be  predicted  as  a  function 
of  the  significant  wave  height  and  the  sediment  fall  velocity.  The 
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other  free  parameters  in  Equation  23  include  the  bed  drag  coefficient 
c^,  the  bedload  efficiency  e^,  and  the  suspended  load  efficiency  e^. 

For  the  present  study  these  variables  were  assumed  to  be  equal  to  0.005, 
0.21,  and  0.025,  respectively.  These  values  were  selected  on  the  basis 
of  an  analysis  of  field  and  laboratory  longshore  transport  data  as 
discussed  by  Bailard  (1981).  Other  similar  values  can  be  selected 
without  qualitatively  changing  the  results. 

Figure  19  shows  a  plot  of  the  predicted  on-offshore  sediment  trans¬ 


port  rate  as  a  function  of  the  significant  wave  height.  The  sediment 
fall  velocity  was  assumed  to  be  equal  to  4  cm/sec,  which  is  appropriate 
for  the  sand  found  at  Torrey  Pines  Beach.  The  bedload-transport  rate  is 
depicted  by  the  dashed  line,  the  suspended  load  transport  rate  by  the 
dotted  line,  and  the  total  load  transport  rate  by  the  solid  line. 

Figure  19  suggests  that,  for  conditions  similar  to  those  at  Torrey  Pines 
Beach  during  the  NSTS  experiment,  sand  is  moved  onshore  when  the  sig¬ 
nificant  wave  height  is  less  than  approximately  90  cm  and  offshore  when 
it  is  greater.  The  maximum  onshore  transport  rate  (0.8  m^m  ^  day 
occurs  when  is  equal  to  approximately  59  cm.  For  waves  with  a  sig¬ 
nificant  wave  height  greater  than  approximately  150  cm,  the  transport 
rate  is  larger  by  a  factor  of  10  and  directed  offshore.  Under  roost 
conditions,  the  predicted  bedload  and  suspended  load  transports  are  both 
in  the  same  direction.  Near  the  null  point  (H^  -  90  cm),  however,  the 
bedload  is  directed  offshore  while  the  suspended  load  is  directed  onshore. 

Qualitatively,  Figure  19  tends  to  confirm  some  aspects  of  observed 
beach  behavior.  During  prolonged  periods  of  small  waves,  a  beach  is 
seen  to  slowly  accrete.  With  the  appearance  of  the  first  large  swell, 
however,  the  beach  can  cut  back  dramatically  within  a  few  days.  Figure  19 
supports  this  observation.  In  addition.  Figure  19  qualitatively  supports 
Short's  (1978)  observations  at  several  Australian  beaches  that  the 
neutral  point  wave  height  separating  accretion  from  erosion  was  equal  to 
120  cm.  Although  the  present  study  suggests  a  neutral  point  wave  height 
of  90  cm,  the  magnitudes  are  similar. 

Because  the  wave  velocity  characteristics  are  independent  of  grain 
size,  the  effect  of  different  sediment  fall  velocities  on  the  total  load 
sediment  transport  rate  may  be  predicted  as  well  (Figure  20).  Increasing 
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the  fall  velocity  is  seen  to  decrease  the  magnitude  of  the  on-offshore 
transport  rate  without  seriously  affecting  the  value  of  the  neutral 
point  wave  height.  As  a  result,  beaches  with  fine  sand  would  be  expected 
to  experience  greater  changes  in  volume  than  equivalent  beaches  with 
coarse  sand. 


DISCUSSION 

The  applicability  of  these  results  to  other  sites  and  to  other  wave 
conditions  is  unknown;  however,  it  is  hypothesized  that  the  results  may 
be  relatively  site  and  time  specific.  Generally  speaking,  for  mono¬ 
chromatic  waves  and  plane  contour  beaches,  the  surfzone  hydrodynamics 
should  be  a  function  of  the  incident  wave  height,  direction,  and  period, 
as  well  as  the  beach  slope.  For  random  waves,  the  shapes  of  the  energy 
and  directional  spectrums  may  also  be  important.  During  the  NSTS  Torrey 
Pines  experiment,  the  waves  were  almost  normally  incident  with  a  near 
uniform  period.  Moreover  the  beach  slope  changed  little  during  the 
month.  The  only  parameter  that  varied  to  any  significant  degree  was  the 
wave  height,  which  varied  by  a  factor  of  3.  As  a  result,  Figures  19  and 
20  cannot  be  directly  extended  to  more  general  conditions.  Neglecting 
for  a  moment  the  effects  of  incident  wave  angle  and  wave  period,  physical 
reasoning  suggests  that  there  should  be  a  family  of  curves  similar  in 
shape  to  those  in  Figures  19  and  20  that  vary  with  beach  slope.  The 
primary  difference  between  different  curves  would  be  the  position  of 
their  neutral  point  wave  height.  For  flat  beaches  typical  of  large  wave 
conditions,  the  neutral  point  wave  height  would  be  larger,  while  for 
steeper  beach  conditions  the  neutral  point  wave  height  would  be  smaller. 
Figures  19  and  20  are  believed  to  represent  transition  conditions. 

Partial  support  for  this  hypothesis  is  provided  in  the  study  by  Aubrey 
(1978),  who  found  that  future  beach  profiles  at  Torrey  Pines  Beach  were 
best  predicted  when  the  existing  beach  profile  shape  and  the  incident 
wave  height  were  known.  Clearly,  more  data  sets  with  a  wider  range  of 
beach  slope  and  wave  conditions  are  needed  to  develop  a  useful  beach 
profile  predictive  capability. 
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Further  discussion  is  also  needed  concerning  possible  errors  in  the 
estimated  surfzone  velocity  moments.  The  odd  moment  quantities  4*^ , 
and  6^,  in  effect,  represent  small  differences  between  large  numbers,  so 
they  are  especially  sensitive  to  small  errors  in  the  measured  current 
meter  data.  Because  of  this  sensitivity,  considerable  care  was  exer¬ 
cised  in  prescreening  the  data.  The  quantity  u*  was  found  to  be  partic¬ 
ularly  sensitive  to  data  errors  and  was  used  as  an  indicator  of  bad 
sections  of  data.  In  spite  of  this  care,  the  data  itself  could  be 
subject  to  an  inherent  bias  due  to  current  meter  inaccuracies. 

One  error  in  particular  may  be  that  the  mean  on-offshore  currents 

were  a  manifestation  of  a  current  meter  rectification  process.*  Some 

evidence,  however,  suggests  that  rectification  may  not  have  been  too 

significant.  First,  if  rectification  were  important  the  mean  offshore 

current  might  be  expected  to  be  a  function  of  u  .  The  latter  was  found 

m 

to  be  relatively  constant  across  the  surfzone,  while  the  mean  offshore 
currents  varied  significantly.  This  tends  not  to  support  a  rectifica¬ 
tion  hypothesis.  Moreover,  Wright  et  al.  (1978)  reported  measuring 
onshore  currents  in  the  upper  part  of  the  water  column  inside  the  surf¬ 
zone  and  offshore  currents  near  the  bottom.  These  measurements  were 
made  with  small  ducted  fan  current  meters  unlike  the  electromagnetic 
current  meters  used  in  the  NSTS  experiments  and,  presumably,  would  not 
be  subject  to  the  same  rectification  characteristics.  Nevertheless, 
until  more  exhaustive  studies  are  done  on  the  response  of  the  electro¬ 
magnetic  current  meters,  the  magnitudes  and  directions  of  the  mean 
on-offshore  currents  must  be  suspect. 

It  is  disappointing  that  the  present  sediment  transport  model  is 
not  more  accurate  at  predicting  the  measured  beach  volume  changes.  The 
reason  for  its  lack  of  accuracy  is  unknown.  The  variability  in  the 
estimated  on-offshore  transport  rates  between  consecutive  hourly  data 
sets  suggests,  however,  that  the  on-offshore  transport  rate  may  vary 
significantly  on  an  hourly  basis  due  to  small  changes  in  the  incident 
wave  field  and  the  tide.  The  latter  changes  the  position  of  the  surf¬ 
zone  relative  to  the  existing  beach  profile,  which  can  significantly 


*From  personal  comnunication  with  R.  T.  Guza  and  D.  G.  Aubrey. 
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alter  the  breaking  wave  characteristics.  As  a  result,  it  may  be  intrin¬ 
sically  difficult  to  test  the  capability  of  an  on-offshore  sediment 
transport  model  to  predict  daily  beach  volume  changes  using  current 
meter  records  of  very  limited  duration. 

In  addition,  some  of  the  measured  beach  volume  changes  seem  to  be 
anomalous  in  light  of  our  present  limited  knowledge.  In  particular,  the 
large  on-offshore  rates  of  sand  movements  during  periods  of  small  waves 
on  November  6,  7,  and  20  are  difficult  to  understand.  It  would  appear 
that  wave  energy  should  have  been  insufficient  to  generate  this  volume 
of  sediment  transport.  One  possible  explanation  may  be  a  temporary 
local  convergence  or  divergence  of  the  longshore  transport  rate.  Another 
explanation  may  be  inaccuracies  in  the  beach  profile  measurements. 


CONCLUSIONS 

Based  on  the  results  of  this  study,  the  following  conclusions  can 
be  made: 


1.  A  thorough  evaluation  of  Bailard's  (1981)  surfzone  sediment 
transport  model  using  daily  beach  volume  measurements  was  not  possible 
due  to  inadequate  surfzone  current  meter  record  lengths  and  to  possible 
inaccuracies  in  the  beach  profile  measurements. 

2.  Predicting  ^each  volume  changes  using  a  highly  Simplified  lora 
of  Bailard's  (1981)  surfzone  sediment  transport  model  appears  promising. 
Analysis  of  additional  surfzone  current  meter  data  sets  having  different 
wave  and  beach  slope  characteristics  is  necessary  before  a  useful  pre¬ 
dictive  capability  can  be  developed. 

3.  Average  surfzone  wave  velocity  moments  were  found  to  be  a 
linear  function  of  significant  wave  height.  Additional  data  sets  are 
needed  to  determine  their  relationships  to  wave  period  and  average  beach 
slope. 
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4.  Except  for  the  mean  on-offshore  current  6  ,  and  the  second 

A  u 

skevmess  parameter  the  normalized  surfzone  distributions  of  the  wave 
velocity  moments  appear  to  exhibit  a  predictable  form  of  behavior. 

5.  Without  analysis  of  additional  data  sets,  the  degree  of  gener¬ 
ality  of  these  results  is  unknown. 


RECOMMENDATIONS 

The  following  research  is  needed  to  develop  a  more  generally  valid 
on-offshore  sediment  transport  model: 

1.  Investigate  pertinent  surfzone  velocity  moments  on  beaches  with 
significantly  different  beach  slopes  and  wave  pe. iods  than  those  encountered 
at  Torrey  Pines  Beach.  The  NSTS  Leadbetter  Beach  data  set  does  provide 

a  beach  with  a  steeper  beach  slope,  so  this  would  be  a  likely  first 
candidate. 

2.  The  significance  of  ignoring  threshold  of  motion  effects  in  the 
present  model  must  be  critically  evaluated. 

3.  A  more  detailed  treatment  of  the  drag  coefficient  may  be 
warranted.  Ideally,  the  coefficient  should  be  estimated  from  longshore 
current  data. 

4.  Sand  tracer  studies  over  a  very  limited  area  (several  meters) 
in  conjunction  with  a  single  current  meter  might  be  beneficial  in  evalu¬ 
ating  potential  sediment  transport  models. 
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Table  1.  Estimated  and  Measured  Beach  Volume  Changes  and  Incident  Wave  Characteristics 


Estimated  as  half  the  2-day  change. 


Date 

(Nov 

1978) 


Table  2.  Average  Surfzone  Velocity  Parameters 
and  Incident  Wave  Characteristics 


f  ^  \ 

(cm  sec  *) 


s 

‘'l 

4* 

2 

0.00126 

0.206 

IQ 

0.00148 

0.197 

mm 

0.00117 

0.234 

IQ 

0.00163 

0.097 

-0.019 

0.00303 

0.106 

-0.152 

0.00341 

0.134 

-0.035 

0.00124 

0.222 

0.326 

0.00194 

0.232 

0.352 

0.00128 

0.223 

0.350 

-0.070 


.550  1.25 

.603  1.20 

.594  1.35 

0.731  1.26 

0.687  1.11 

0.574  0.940 

0.587  1.27 

0.574  1.24 

0.574  1.28 


Derme  Pirimetcn 


Predicted  normalized  equilibrium  beach  slope,  tan  p/ak^,  as  a  function  of  the  normalized 
depth  k„h  and  the  wave  amplitude  parameter  aa/W  (from  Bailard,  1981 ). 


2-dircctional  current  meters 


Figure  4.  Relative  current  meter  and  pressure  sensor  positions  for  the  NSTS  Torrey  Pines 
.  experiment  (adapted  from  Seymour  and  King,  1982). 


estimated  bedload  and  suspended  load 


Figure  6.  Average  surfzone  values  of  the  skewness  parameter  versus  significant  wave 
heights  computed  from  single  64-minute  current  meter  records.  The  straight 
line  is  a  least-squares  fit  to  the  data. 


IS  significant  wave 
It  line  is  a  least- 


Figure  8.  Average  surfzone  values  of  the  normalized  mean  onshore  current  6^  versus  significant 
wave  height  computed  from  single  64-minute  records.  The  straight  line  is  a  least- 
squares  fit  to  the  data. 


Average  surfzone  values  of  the  orbital  velocity  magnitude  Ujj,  versus  significant 
wave  height  computed  from  single  64-minute  records.  The  straight  line  is  a 
least-squares  fit  to  the  data. 
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Figure  12.  Average  surfzone  values  of  the  standard  deviation  ^2^d.  versus  significant 
wave  height  computed  from  single  64-minute  records. 
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te  that  a  positive  normalized  current  is  directed  offshore. 


Normalized  Surfzonc  Position, 


Significant  Wave  Height,  Hj  (cm) 

Figure  20.  Predicted  total  load  on-offshore  sediment  transport  rate  as  function  of  significant  wave 
height  and  sediment  fall  velocity.  Greater  rates  of  transport  are  found  for  smaller 
sized  sediments. 
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BUREAU  OF  RECLAMATION  Code  L‘il2  (C  Se(ander)  Denver  CO 

CINCLANT  CIV  ENGR  SUPP  PLANS  OFFR  NORFOLK.  VA 

CINCPAC  Fac  Engrng  Div  (J44)  Makalapa.  HI 

CNM  Code  0.2462,  Washington  DC;  Code  043  Washington  DC;  Code  MAT-04,  Washington.  DC;  MAT-0718. 
Washington,  DC;  NMAT  -  044.  Washington  DC 

CNO  Code  NOP-964,  Washington  DC;  Code  OP  323,  Washington  DC;  Code  OP  40.S.  Washington  DC;  Code  OP 
405.  Washington.  DC;  Code  OP  414.  Washington  DC;  Code  OP  97  Washington  DC;  Code  OP  97 
Washington,  DC:  Code  OP  987  Washington  DC;  Code  OP323  Washington  DC;  Code  OPNAV  09B24  (H); 
Code  OPNAV  22.  Wash  DC;  Code  OPNAV  23.  Wash  DC;  OP-41  IF.  Wash  DC;  OP987J.  Washington.  DC 
COMCBPAC  Operations  Off.  Makalapa  HI 
COMNAVBEACHPHIBREFTRAGRU  ONE  San  Diego  CA 
COMNAVMARIANAS  Code  N4,  Guam 
COMNAVSURFLANT  Norfolk.  VA 
COMOCEANSYSPAC  SCE.  Pearl  Harbor  HI 
COMRNCF  Nicholson,  Tampa.  FL;  Nicholson.  Tampa.  FL 
COMSUBDEVGRUONE  Operations  Offr,  San  Diego,  CA 
NAVSURFPAC  Code  N-4.  Coronado 

COMOPTEVFOR  CMDR.  Norfolk.  VA;  Code  705.  San  Diego,  CA 

DEFFUELSUPPCEN  DFSC-OWE  (Term  Engrng)  Alexandria,  VA;  DFSC-OWE.  Alexandria  VA 

DLSIE  Army  Logistics  Mgt  Center,  Fort  Lee,  VA 

DNA  STTL,  Washington  DC 

DTIC  Defense  Technical  Info  Ctr'Alexandria,  VA 

DTNSRDC  Anna  Lab  (Code  1175)  Annapolis  MD;  Anna  Lab  (Code  119)  Annapolis  MD;  Anna  Lab  (Code 
1568)  Annapolis  MD 
FMFLANT  CEC  Offr.  Norfolk  VA 
FMFPAC  CG(FEO)  Camp  Smith,  HI 

GSA  Assist  Comm  Des  &  Cnst  (FAI.A)  D  R  Dibner  Washington.  DC 
LIBRARY  OF  CONGRESS  Washington,  DC  (Sciences  &  Tech  Div) 

MARINE  CORPS  BASE  Isl  For  Scrv  Supp  Gru  (CSS-5)  Camp  Pendleton  CA;  PWO,  Camp  Pendleton  CA; 
PWO.  Camp  S.  D.  Butler,  Kawasaki  Japan 

MCAS  Facil.  Engr.  Div.  Cherry  Point  NC;  Facs  Maint  Dept  -  Operations  Div,  Cherry  Point 

MCDEC  M&L  Div  Ouantico  VA;  NSAP  REP.  Ouantico  VA 

MCRD  SCE.  San  Diego  CA 

MILITARY  SEALIFT  COMMAND  Washington  DC 

NAF  PWO.  Atsugi  Japan 

NARF  Equipment  Engineering  Division  (Code  61000),  Pensacola.  FL 

NAS  Dir  of  Engrng,  PWD.  Corpus  Christi.  TX;  Dir,  Util.  Div.,  Bermuda;  Lead.  Chief.  Petty  Offr.  PW/Self 
Help  Div.  Beevilic  TX;  PW  (J.  Maguire).  Corpus  Christi  TX;  PWD  -  Engr  Div.  Oak  Harbor,  WA;  PWD 
Maint.  Div..  New  Orleans.  Belle  Chasse  LA;  PWD,  Code  I82IH  (Pfankuch)  Miramar.  SD  CA;  PWO  Belle 
Chasse,  LA;  PWO  Key  West  FL;  PWO.  Glenview  IL 
NAVFACENGCOM  CONTRACTS  ROlLC  Key  West  FL 
NAS  SCE  Nc  'k.  VA 

NATL  BUREAU  OF  STANDARDS  Kovacs,  Washington.  D.C,;  R  Chung  Washington.  DC 
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NATl.  RESEARCH  COL’NCIL  Naval  Studies  Board,  Washington  D< 

NAVACT  PWC).  London  L'K 
NAVAEROSPREGMEDCEN  SCE,  Pcnsavola  FI 
NAVAIRDEVCEN  Code  8LV  W'arminster  PA 

NAVCHAPGRC  CO  Williamsburg  \  .A,  0|veralions  Olluct.  (  ode  (o  Willi.irisbuie  \A 
NA\’COASTS^’SCE.N  CO,  Panama  ('itv  FL,  Code  "'1^  IJ  Ouirk)  Panama  (iiv.  FI,,  (ode  715  (J  Mittleman) 
Panama  City,  FL:  Code  ’’77  IC  B  Koesyl  Panama  Cits  FI  :  Library  Panama  (ns.  FL.  PWO  Panama  City.  FL 
NAVCOMMAREA.V1STRSIA  S(  E  Cnit  I  Naples  Italy;  SCE  Wahiassa  HI 
NAVCOMMSTA  Code  -Ull  Nea  Makri.  Greece,  PWO.  1  xmouih.  Australia 
NAVCONSTRACEN  Curriculum  Instr  Stds  Offr.  Gulfport  MS 
NAVEDTRAPRODEVCEN  Technical  Librars .  Pens., cola.  FI 

NAVELEXSVSCCFM  Code  PME  i:4.M,  Washington.  DC  ,  PML  i:4-(si:.  Wash  IK 
NA\  EODTECHCEN  Code  «I5.  Indian  Head  MD 
NANFAC  PWO.  Centersillu  Bch,  Fcrndalc  CA 

N.AVF.ACENGC'OM  Code  (I.4T  (Essoglou)  Alexandria.  V’A.  Code  (14.4  Alexandria.  \'.A.  (  ode  1144  Alexandria. 

VA.  Code  t)44BI  Alexandria.  VA;  Code  (14.51  (P  W  Brewer)  Alexandria.  Va,  Code  045.4  iD  Potter) 
Alexandria.  VA,  Code  l)4.5.4C,  Alexandria.  VA;  Code  (MAI  Alexandria.  V.A;  Code  Ob.  Alexandria  \'A; 

Code  0‘)M54.  Technical  Library,  Alexandria.  VA;  Code  111.4.  Alexandria.  VA 
NAVFACENGCOM  CONTRACTS  ROICC.  4  ap:  ROICC  Code  495  Portsmouth  VA 
NAVFACENGC'OM  -  CHES  DIV.  Code  405  Wash.  DC;  Code  407  (D  Scheesele)  Washington.  DC;  Code 
FTO-IC  Washington  DC;  FPO-I  W'ashmgton,  DC';  FPO-IEA5  Washington  DC;  Library.  Washington,  D  C 
NAVFACENGCOM  -  LANT  DIV.  Code  405  Civil  Engr  BR  Norfolli  VA;  Eur  BR  Deputy  Dir,  Naples  Italy; 
Ltbrary.  Norfolk,  VA;  RDT&ELO  102A,  Norfolk,  VA 

NAVFACENGCOM  -  NORTH  DIV  (Boretsky)  Philadelphia.  PA.  CO.  Code  09P  Philadelphia  PA;  Code  1028. 

RDT&ELO,  Philadelphia  PA;  Library.  Philadelphia.  PA;  ROICC.  Contracts.  Crane  IN 
NAVFACENGCOM  -  PAC  DIV.  CODE  09P  PEARL  HARBOR  HI.  Code  2011  Pearl  Harbor,  HI;  Code  402. 

RDT&E.  Pearl  Harbor  HI;  Commander.  Pearl  Harbor.  HI;  Librars.  Pearl  Harbor.  HI 
NAVFACENGCOM  -  SOUTH  DIV.  Code  90,  RDT&ELO.  Charleslon  SC;  Librars,  Charleston,  SC 
NAVFACENGCOM  -  WEST  DIV.  Code  04B  San  Bruno.  CA;  Library.  San  Bruno.  CA;  O9P'20  San  Bruno. 

CA;  RDT&ELO  Code  2011  San  Bruno.  CA 

NAVFACENGCOM  CONTRACTS  Eng  Div  dir.  Southwest  Pac.  Manila.  PI;  OICC.  Southwest  Pac.  Manila.  PI; 

OICC/ROICC.  Balboa  Panama  Canal;  ROICC.  Keflavik.  Iceland:  ROICC.  NAS.  Corpus  Christi.  TX 
NAVFORCARIB  Commander  (N42),  Puerto  Rico 
NAVMAG  SCE.  Subic  Bay,  R.P. 

NAVOCEANO  Library  Bay  St.  Louis.  MS 

NAVOCEANSYSCEN  Code  4473  Bayside  Library,  San  Diego.  CA;  Code  4473B  (Tech  Lib)  San  Diego.  CA; 

Code  09  (Talkington),  San  Diego.  CA;  Code  5204  (J.  Stachiw).  San  Diego.  CA;  Code  5214  (H.  Wheeler). 

San  Diego  CA;  Code  5221  (R. Jones)  San  Diego  Ca;  Code  5311  (Bachman)  San  Diego.  CA;  Hawaii  Lab  (R 
Yumori)  Kailua.  HI;  Hi  Lab  Tech  Lib  Kailua  HI 

NAVPGSCOL  C.  Morers  Monterey  CA;  Code  61WL  (O.  Wilson)  Monterey  CA;  E.  Thornton.  Monterey  CA 
N.AVPHIBASE  CO,  ACB  2  Norfolk.  VA;  COMNAVBEACHGRU  TWO  Norfolk  VA;  Code  S3T.  Norfolk  VA; 

Dir.  Amphib.  Warfare  Brd  Staff,  Norfolk,  VA;  Harbor  Clearance  Unit  Two.  Little  Creek.  VA 
NAVFACENGCOM  CONTRACTS  OICC/ROICC,  Norfolk,  VA 
NAVPHIBASE  SCE  Coronado,  SD.CA 
NAVREGMEDCEN  SCE;  SCE,  Guam 
NAVSCOLCECOFF  C35  Port  Hueneme.  CA 
NAVSeSOL  PWO,  Athens  GA 

NAVSEASYSCOM  Code  0325,  Program  Mgr,  Washington,  DC;  Code  SEA  OOC  Washington,  DC;  SEA  04E  (L 
Kess)  Washington,  DC 
NAVSECGRUACT  PWO,  Adak  AK 

NAVSHIPYD  Bremerton,  WA  (Carr  Inlet  Acoustic  Range);  Code  202.5  (Library)  Puget  Sound.  Bremerton 
WA;  Code  380,  Portsmouth,  VA;  Code  440  Portsmouth  NH;  Code  440,  Puget  Sound,  Bremerton  WA;  Tech 
Library,  Vallejo,  CA 

NAVSTA  Dir  Engr  Div,  PWD,  Mayport  FL;  PWD  (LTJG.P.M.  Motolenich).  Puerto  Rico;  PWO,  Mayport  FL; 

Utilities  Engr  Off.  Rota  Spain 
NAVTECHTRACEN  SCE,  Pensacola  FL 

NAVWPNCEN  Code  2636  China  Uke;  Code  3803  China  Uke,  CA 
NAVWPNSTA  PW  Office  Yorktown,  VA 

NAVWPNSTA  PWD  -  Maim.  Control  Div.,  Concord,  CA;  PWO.  Charleston,  SC;  PWO.  Seal  Beach  CA 
NAVWPNSUPPCEN  Code  09  Crane  IN 
NCTC  Const.  Elec.  School,  Port  Hueneme,  CA 

NCBC  Code  10  Davisville,  Rl;  Code  15,  Port  Hueneme  CA;  Code  155.  Port  Hueneme  CA;  Code  156,  Port 
Hueneme,  CA;  Code  1571,  Port  Hueneme,  CA 
NCBU  411  OIC,  Norfolk  VA 
NCR  20,  Commander;  30th  Det,  OIC,  Diego  Garcia  I 
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NMCB  74.  CO;  FI\'E.  Operations  Dept 

NOAA  (Dr  T  Me  Guinness)  RiK'kville.  MD;  Library  RtK'kvillc.  MD 
NORDA  Code  41tl  Bay  St  Louis.  MS.  Code  441)  (Ocean  Rsch  Off)  Bay  St.  Louis  MS 
NRL  Code  .S8))l)  Washington.  DC.  Code  .‘'S4.f  (F  Rosenthal)  Washington.  DC;  Code  8441  (R.A.  Skop), 
Washington  DL' 

NROTC  LW  Stephenson.  L'C  Berkeley.  CA 

NSD  SCE.  Subic  Bay.  R  P 

NTC  OICC,  CBU-4()1.  Great  Lakes  IL 

NL'CLEAR  REGULATORY  COMMISSION  T  C.  Johnson,  Washington.  DC 

NUSC  Code  l.Jl  New  London.  CT;  Code  332.  B-80  (J  Wilcox)  Newt  London.  CT;  Code  EA123  (R.S.  Munn). 

New  I^ndon  CF;  Code  TA1.31  (G.  De  la  Cruz).  New  London  CT 
OFFICE  SECREI  ARY  OF  DEFENSE  ASD  (MRA&L)  Code  CSS/CC  Washington.  DC 
ONR  Central  Regional  Office,  Boston.  MA;  Code  481,  Bay  St.  Louis,  MS;  Code  485  (Silva)  Arlington,  VA; 
Code  700F  Arlington  VA 

PERRY  OCEAN  ENG  R  Pellen.  Riviera  Beach.  FL 

PHIBCB  1  P&E.  San  Diego.  CA;  I,  CO  San  Diego,  CA.  I.  CSWC  D  Wellington.  San  Diego,  CA 
PMTC  EOD  Mobile  Unit.  Point  Mugu.  CA;  Pat  Counsel.  Point  Mugu  CA 

PWC  CO  Norfolk,  VA;  CO,  (Code  10).  Oakland,  CA;  CO.  Great  Lakes  IL;  CO,  Pearl  Harbor  HI;  Code  10. 
Great  Lakes.  IL;  Code  120,  Oakland  CA;  Code  120C.  (Library)  San  Diego.  CA;  Code  128,  Guam;  Code 
154.  Great  Lakes.  IL,  Code  400,  Great  Lakes.  IL;  Code  400,  Pearl  Harbor.  HI;  Code  400,  San  Diego.  CA; 
Code  420.  Great  Lakes.  IL;  Code  420,  Oakland.  CA;  Code  424.  Norfolk.  VA;  Code  700,  Great  Lakes.  IL; 
Code  7(X).  San  Diego.  CA,  Commanding  Officer.  Great  Lakes.  IL;  Library,  Guam;  Library,  Norfolk.  VA; 
Library.  Oakland.  CA;  Library.  Subic  Bay.  R.P..  Library.  Pearl  Harbor.  HI 
SUPANX  PWO,  Williamsburg  V  A 
UCT  ONE  OIC,  Norfolk.  VA 
ucr  rwo  OIC,  Port  Hueneme  CA 

U  S.  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

US  DEPT  OF  INTERIOR  Bur  of  Land  Mgmnt  Code  583,  Washington  DC 

US  GEOLOGICAL  SURVEY  Off  Marine  Geology.  Piteleki.  Reston  VA 

US  NATIONAL  MARINE  FISHERIES  SERVICE  Highlands  NY  (Sandy  Hook  Ub-Library) 

US  NAVAL  FORCES  Korea  (ENJ-P&O) 

USCG  (G-MP-3  USP/82)  Washington  Dc;  (Smith).  Washington.  DC 
USCG  R&D  CENTER  CO  Groton.  CT;  D.  Motherway,  Groton  <T 
L'SDA  Forest  Service  Reg  .3  (R.  Brown)  Albuquerque,  NM 

USNA  ENGRNG  Div.  PWD,  Annapolis  MD;  USNA/'Sys  Eng  Dept.  Annapolis,  MD 
WATER  &  POWER  RESOURCES  SERVICE  (Smoak)  Denver.  CO 
BERKELEY  PW  Engr  Div,  Harrison,  Berkeley,  CA 

CALIF.  DEPT  OF  NAVIGATION  &  OCEAN  DEV.  Sacramento.  CA  (G.  Armstrong) 

CALIF  MARITIME  ACADEMY  Vallejo.  CA  (Library) 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY  Pasadena  CA  (Keck  Ref.  Rm) 

CALIFORNIA  STATE  UNIVERSITY  (Yen)  Long  Beach,  CA;  LONG  BEACH.  CA  (CHELAPATl) 
CATHOLIC  UNIV.  Mech  Engr  Dept.  Prof.  Niedzwecki,  Wash.,  DC 
COLORADO  STATE  UNIV  .  FOOTHILL  CAMPUS  Fort  Collins  (Nelson) 

CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.);  Ithaca.  NY  (Civil  &  Environ.  Engr) 

DAMES  &  MOORE  LIBRARY  LOS  ANGELES,  CA 

DUKE  UNIV  MEDICAL  CENTER  B.  Muga,  Durham  NC;  DURHAM.  NC  (VESIC) 

FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton.  FL  (McAllister) 

HARVARD  UNIV.  Dept,  of  Architecture.  Dr.  Kim.  Cambridge.  MA 
GEORGIA  INSTITUTE  OF  TECHNOLOGY  Atlanta  GA  (B,  Mazanti) 

INSTITUTE  OF  MARINE  SCIENCES  Morehead  City  NC  (Director) 

IOW.\  STATE  UNIVERSITY  Ames  I A  (CE  Dept.  Handy) 

woe  OS  HOLE  OCEANOGRAPHIC  INST.  Woods  Hole  MA  (Winget) 

LEHIGH  UNIVERSITY  BETHLEHEM,  PA  (MARINE  GEOTECHNICAL  LAB.,  RICHARDS);  Bethlehem 
PA  (Fritz  Engr.  Lab  No,  13,  Beedle);  Bethlehem  PA  (Linderman  Lib.  No. .30.  Flecksteiner) 

MAINE  MARITIME  ACADEMY  CASTINE.  ME  (LIBRARY) 

MICniGAN  TECHNOLOGICAL  UNIVERSITY  Houghton,  MI  (Haas) 

MIT  Cambridge  MA;  Cambridge  MA  (Rm  10-500.  Tech.  Reports.  Engr,  Lib.);  Cambridge  MA  (Whitman); 
Cambridge,  MA  (Harleman) 

NAIT  ACADEMY  OF  ENG.  ALEXANDRIA,  VA  (SEARLE,  JR.) 

NATURAL  ENERGY  LAB  Library.  Honolulu,  HI 

NEW  MEXICO  SOLAR  ENERGY  INST.  Dr.  Zwibel  Las  Cruces  NM 

OREGON  STATE  UNIVERSITY  (CE  Dept  Grace)  Corvallis.  OR;  CORVALLIS,  OR  (CE  DEPT.  BELL); 
Corvalis  OR  (School  of  Oceanography) 

PENNSYLVANIA  STATE  UNIVERSITY  STATE  COLLEGE.  PA  (SNYDER);  State  College  PA  (Applied 
Rsch  Lab);  UNIVERSITY  PARK.  PA  (GOTOLSKI) 
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PURDUE  UNIVERSITY  Lafayette  IN  (Leonards):  Lafayette.  IN  (Altschaeffl);  Lafayette.  IN  (CE  Engr.  Lib) 

SAN  DIEGO  STATE  UNIV.  1.  Noorany  San  Diego.  CA 

SCRIPPS  INSTITUTE  OF  OCEANOGRAPHY  LA  JOLLA.  CA  (ADAMS) 

SEATTLE  U  Prof  Schwaegler  Seattle  WA 

SOUTHWEST  RSCH  INST  King.  San  Antonio.  TX;  R.  DeHart.  San  Antonio  TX 
STATE  UNIV.  OF  NEW  YORK  Buffalo,  NY;  Fort  Schuyler,  NY  (Longobardi) 

TEXAS  A&M  UNIVERSITY  College  Station  TX  (CE  Dept.  Herbich) 

TEXAS  TECH  UNIVERSITY  Dept  of  IE  (Prof.  Ayoub),  Lubbock  TX 

UNIVERSITY  OF  ALASKA  Doc  Collections  Fairbanks,  AK;  Marine  Science  Inst.  College.  AK 
UNIVERSITY  OF  CALIFORNIA  A-031  (Storms)  La  Jolla.  CA;  BERKELEY,  CA  (CE  DEPT,  GERWICK); 
BERKELEY,  CA  (CE  DEPT,  MITCHELL);  Berkeley  CA  (Dept  of  Naval  Arch  );  Berkeley  CA  (E. 
Pearson);  La  Jolla  CA  (Acq.  Dept,  Lib.  C-075A);  M.  Duncan,  Berkeley  CA 
UNIVERSITY  OF  CONNECTICUT  Groton  CT  (Inst  Marine  Sci,  Library) 

UNIVERSITY  OF  DELAWARE  Newark,  DE  (Dept  of  Civil  Engineering.  Chesson) 

UNIVERSITY  OF  HAWAII  HONOLULU,  HI  (SCIENCE  AND  TECH.  DIV.) 

UNIVERSITY  OF  ILLINOIS  (Hall)  Urbana.  IL,  Metz  Ref  Rm.  Urbana  IL;  URBANA.  IL  (DAVISSON); 
URBANA.  IL  (LIBRARY) 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus),  ME  Dept,  Amherst,  .MA 
UNIVERSITY  OF  MICHIGAN  Ann  Arbor  MI  (Richart) 

UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln.  NE  (Ross  Ice  Shelf  Pro)  ) 

UNIVERSITY  OF  NEW  HAMPSHIRE  DURHAM.  NH  (LAVOIE) 

UNIVERSITY  OF  SO.  CALIFORNIA  Univ  So.  Calif 

UNIVERSITY  OF  TEXAS  Inst.  Marine  Sci  (Library).  Port  Arkansas  TX 

UNIVERSITY  OF  TEXAS  AT  AUSTIN  Austin  TX  (R.  Olson) 

UNIVERSITY  OF  WASHINGTON  Seattle  WA  <M  Sherif);  SEATTLE.  WA  (APPLIED  PHYSICS  LAB); 
SEATTLE,  WA  (OCEAN  ENG  RSCH  LAB,  GRAY);  SEATTLE.  WA  (PACIFIC  MARINE  ENVIRON 
LAB..  HALPERN);  Seattle  WA  (E.  Linger);  Seattle.  WA  Transportation.  Construction  &  Geom.  Div 
UNIVERSITY  OF  WISCONSIN  Milwaukee  W1  (Ctr  of  Great  Lakes  Studies) 

VIRGINIA  INST.  OF  MARINE  SCI  Gloucester  Point  VA  (Library) 

ALFRED  A.  YEE  &  ASSOC.  Librarian,  Honolulu.  HI 
AMETEK  Offshore  Res.  &  Engr  Div 
ANALYTICAL  TECH  Lawrence.  Port  Hueneme.  CA 
ATLANTIC  RICHFIELD  CO.  DALLAS,  TX  (SMITH) 

BECHTEL  CORP.  SAN  FRANCISCO,  CA  (PHELPS) 

BETHLEHEM  STEEL  CO.  Dismuke.  Bethelehem,  PA 
BRAND  INDUS  SERV  INC,  J.  Buehler,  Hacienda  Heights  CA 
BRITISH  EMBASSY  M  A  Wilkins  (Sci  &  Tech  Dept)  Washington.  DC 
BROWN  &  ROOT  Houston  TX  (D.  Ward) 

CHEVRON  OIL  FIELD  RESEARCH  CO  LA  HABRA,  CA  (BROOKS) 

CICB  O  Rourke 

COLUMBIA  GULF  TRANSMISSION  CO.  HOUSTON,  TX  (ENG.  LIB.) 

CONCRETE  TECHNOLOGY  CORP.  TACOMA,  WA  (ANDERSON) 

DILLINGHAM  PRECAST  F.  McHale,  Honolulu  HI 
DIXIE  DIVING  CENTER  Decatur.  GA 
DRAVO  CORP  Pittsburgh  PA  (Wright) 

EVALUATION  ASSOC.  INC  KING  OF  PRUSSIA,  PA  (FEDELE) 

EXXON  PRODUCTION  RESEARCH  CO  Houston,  TX  (Chao) 

FORD.  BACON  &  DAVIS.  INC.  New  York  (Library) 

FURGO  INC.  Library,  Houston,  TX 

GEOTECHNICAL  ENGINEERS  INC.  Winchester,  MA  (Paulding) 

HALEY  &  ALDRICH,  INC.  Cambridge  MA  (Aldrich.  Jr.) 

JAMES  CO,  R.  Girdley,  Orlando  FL 

LAMONT-DOHERTY  GEOLOGICAL  OBSERVATORY  (Selwyn)  Palisades.  NY;  Palisades  NY  (McCoy) 
LOCKHEED  MISSILES  &  SPACE  CO.  INC.  L.  Trimble,  Sunnyvale  CA 
MARATHON  OIL  CO  Houston  TX 

MARINE  CONCRETE  STRUCTURES  INC.  MEFAIRIE,  LA  (INGRAHAM) 

MC  CLELLAND  ENGINEERS  INC  Corp  Library  Houston.  TX 
MOBIL  PIPE  LINE  CO.  DALLAS.  TX  MGR  OF  ENGR  (NOACK) 

MOFFATT  Si  NICHOL  ENGINEERS  (R.  Palmer)  Long  Beach,  CA 

NEWPORT  NEWS  SHIPBLDG  &  DRYDOCK  CO.  Newport  News  VA  (Tech,  Lib.) 

OPPENHEIM  Los  Angeles.  CA 

PACinC  MARINE  TECHNOLOGY  (M.  Wagner)  Duvall.  WA 

PORTLAND  CEMENT  ASSOC.  SKOKIE,  IL  (CORLEY;  Skokie  IL  (Rsch  &  Dev  Ub,  Lib.) 

R  J  BROWN  ASSOC  (McKeehan).  Houston.  TX 

RAYMOND  INTERNATIONAL  INC.  E  Colie  Soil  Tech  Dept,  Pennsauken,  NJ;  J.  Welsh  Soiltech  Dept. 
Pennsauken,  NJ 


SANDIA  LABORATORIES  Library  Div.,  Livermore  CA;  Seabed  Progren  Div  4536  (D.  Talbert)  Albuquerque 
NM 

SCHUPACK  ASSOC  SO.  NORWALK,  CT  (SCHUPACK) 

SEATECH  CORP,  MIAMI,  FL  (PERONI) 

SHANNON  &  WILLSON  INC.  Librarian  Seattle,  WA 

SHELL  DEVELOPMENT  CO.  Houston  TX  (C.  Sellars  Jr  );  Houston  TX  (E.  Doyle) 

SHELL  OIL  CO.  HOUSTON,  TX  (MARSHALL);  Houston  TX  (R.  de  Castongrene);  I.  Boaz,  Houston  TX 
TIDEWATER  CONSTR.  CO  Norfolk  VA  (Fowler) 

WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib.  Bryan) 

WESTINTRUCORP  Egerton,  Oxnard,  CA 

WM  CLAPP  LABS  -  BATTELLE  DUXBURY,  MA  (LIBRARY) 

WOODWARD-CLYDE  CONSULTANTS  (A.  Harrigan)  San  Francisco;  PLYMOUTH  MEETING  PA  (CROSS, 
III) 

AL  SMOOTS  Los  Angeles,  CA 

ANTON  TEDESKO  Bronxville  NY 

BARA,  JOHN  P.  Lakewood,  CO 

BRAHTZ  La  Jolla,  CA 

BULLOCK  La  Canada 

F.  HEUZE  Alamo,  CA 

GERWICK,  BEN  C.  JR  San  Francisco,  CA 

LAYTON  Redmond,  WA 

CAPT  MURPHY  Sunnyvale,  CA 

PAULI  Silver  Spring,  MD 

R.F.  BESIER  Old  Saybrook  CT 

BROWN  &  CALDWELL  Saunders,  E.M./Oakland,  CA 

WM  TALBOT  Orange  CA 
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